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The transmission spectrum of a high-finesse optical cavity containing an arbitrary number of 
trapped atoms is presented. We take spatial and motional effects into account and show that in the 
limit of strong coupling, the important spectral features can be determined for an arbitrary number 
of atoms, N. We also show that these results have important ramifications in limiting our ability 
to determine the number of atoms in the cavity. 

I. INTRODUCTION 

Cavity quantum electrodynamics (CQED) in the strong coupling regime holds great interest for experimentalists 
and theorists for many reasons From an applied perspective, CQED provides precise tools for the fabrication 

of devices which generate useful output states of light, as exemplified by the single-photon source 0,0, @|i the N- 
photon source Q, and the optical phase gate y|. Conversely, CQED effects transform the high- finesse cavity into a 
sensitive optical detector of objects which are in the cavity field. Viewed simply, standard optical microscopy is made 
more sensitive by having a probe beam pass through the sample multiple times, and by efficiently collecting scattered 
light. In the weak-coupling regime, this has allowed for nanometer- resolution measurements of the positions of a 
trapped ion [^.llOfl. In the strong-coupling regime, the presence and position of single atoms can be detected with high 
sensitivity by monitoring the transmission IT^ , phase shift [l3l | , or spatial mode |14| of probe light sent through 
the cavity. 

In this paper, we consider using strong-coupling CQED effects to precisely count the number of atoms trapped 
inside a high-finesse optical microcavity. The principle for such detection is straightforward: the presence of atoms 
in the cavity field splits and shifts the cavity transmission resonance. A precise iV-atom counter could be used to 
prepare the atoms-cavity system for generation of optical Fock states of large photon number Q , or to study ultra-cold 
gaseous atomic s y ste ms g ases in which atom number fluctuations are important, such as number-squeezed |l5j and 
spin-squeezed \lQL Il7t [i"8| systems. 

A crucial issue to address in considering such a CQED device is the role of the spatial distribution of atoms and their 
motion in the cavity field. An N-atom counter (or any CQED device) would be understood trivially if the N atoms 
to be counted were held at known, fixed positions in the cavit y fie ld. This is a central motivation for the integration 
of CQED with extremely strong traps for neutral atoms [Hi l20j or ions 

UGH- The Tavis-Cummings model H2, 
which applies to this case, predicts that the transmission spectrum of a cavity containing N identically-coupled (with 
strength g), resonant atoms will be shifted from the empty cavity resonance by a frequency g^/N at low light levels. 
Atoms in a cavity can then be counted by measuring the frequency shift of the maximum cavity transmission, and 
distinguishing the transmission spectrum of N atoms from that of N + 1 atoms in the cavity. However, to assess the 
potential for precise CQED-aided probing of a many-body atomic system, we consider here the possibility that atoms 
are confined at length scales comparable to or indeed larger than the optical wavelength. 

In this paper, we characterize the influence of cavity mode spatial dependence and atomic motion on the transmission 
spectrum for an arbitrary number of atoms. The impact of atomic motion on CQED has been addressed theoretically 
in previous work p3 . 123 . |24| , although attention has focused primarily on the simpler problem of a single atom in 
the cavity field. We show that when spatial dependence is included, the intrinsic limits on atom counting change 
significantly. The organization of this paper as follows. In Sec. [n] we introduce the system Hamiltonian, define 
our notation, and derive an explicit expression for the intrinsic transmission function. In Sec. IIIII we introduce the 
method of moments, and use this method to calculate the shape of the intrinsic transmission function. Conclusions 
and implications for atom counting are presented in Sec. IIVI 

II. TRANSMISSION 

Let us consider the Hamiltonian for N identical two-level atoms in a harmonic potential inside an optical cavity 
which admits a single standing wave mode of light. We consider atomic motion and the spatial variation of the cavity 
mode only along the cavity axis, assuming that the atoms are confined tightly with respect to the cavity mode waist 
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in the other two dimensions. The Hamiltonian for this system is 

N 

H = hcj c a1a + J2 huJ a\e l ) (ei\ + H + V (1) 

i=l 

where oj c is the frequency of the cavity mode and a(a^) is the annihilation (creation) operator for the cavity field. 
The motional Hamiltonian Hq = X^-^0,i is a sum over single-atom Hamiltonians Hg t i = pf/2m + mujQxf/2 where 
?7i is the atomic mass and u>o the harmonic trap frequency. The atomic ground and excited internal states, \g) and 
|e), respectively, are separated by energy Tioj a . The dipole interaction with the light field V — J2i Vi is a sum over 
interactions with the dipole moment of each atom Vi = fig cos(kxi) (jej (gi\ a + \gi) (e^ | a)\ where g is the vacuum- 
Rabi splitting, which depends on the atomic dipole moment and the volume of the cavity mode. In this paper we 
assume the cavity mode frequency to be in exact resonance with the atomic transition frequency, uj c — uj a . 

Since the Hamiltonian (Eq. ^| commutes with the total excitation operator, tit = a)a + \a) (ej|, the eigen- 
spectrum of H breaks up into manifolds labelled by their total excitation number. In this work, we are concerned 
with excitation spectra of the atoms-cavity system at the limit of low light intensity, and we therefore restrict our 
treatment to the lowest two manifolds, with = {0, 1}. 

In particular, we consider the excitation spectra from the ground state (motional and internal) of the atoms-cavity 
system. This state |\&o) is given simply as a product of motional and internal states, |^o) = |$i)®|0 c ; gi, g%, ■ ■ . pjv)- In 
the uncoupled internal state notation, the C symbol indicates there are zero photons in the cavity, and the gi symbol 
indicates that atom i is in the ground state. The motional state = TiiLi l<?^o (^i )) is a product of single-atom 
ground states of the harmonic trap. 

Let us calculate the low-light intensity transmission spectrum of the cavity. We assume that the system is pumped by 
a near-resonant linearly coupled driving field such that the cavity excitation Hamiltonian is Hj = E (a^e~ laJt + ae laJt ) 
where E is the product of the external driving electric field strength and the transmissivity of the input cavity mirror, 
and u> is the driving frequency. To determine the cavity transmission spectrum, we determine the excitation rate 
to atoms-cavity states in the tit = 1 manifold from the initial ground state. The atoms-cavity eigenstates decay 
cither by cavity emission, with the transmitted optical power proportional to k (N c ) where k is the cavity decay rate 
and N c = a^a is the intracavity photon number operator, or by other processes (spontaneous emission, losses at the 
mirrors, etc.) at the phenomenological rate constant 7. Neglecting the width of the transmission spectrum caused by 
cavity and atomic decay (k, 7 — > 0), we use Fermi's Golden Rule to obtain the transmission spectrum I(u>): 

/Hex Y, |<*i|a + l*o>r% 3 --wo-o;)= £ |<%|*j)| 2 Sfa - u> - w), (2) 

j,n T = l j,n T = l 

where = a' \^o)- In the summation over all atoms-cavity eigenstates, we make the simplification that only states 
with tit — 1 need be included since only these states are coupled to the ground state by a single excitation. To 
simplify notation, we make this implicit assumption throughout the remainder of this paper. We denote I(ui) as the 
"intrinsic transmission spectrum" . In the limit of k, 7 — > this is composed of delta functions in frequency, while an 
experimentally observed transmission spectrum would be convolved by non-zero linewidths. 

To proceed further, we introduce the basis states {|0) ; \i)} which span the space of internal states in the — 1 
manifold. The state |0) = \l c \ 9o,9i, ■ ■ ■ 9n) has one cavity photon and all atoms in their ground state. The state 
\i) = |0 C ; 5oi3ii ■ 9n) is the state in which the cavity field is empty, while a single atom (atom i) is in the 

excited state. Restricted to the ut — 1 manifold, the Hamiltonian (Eq. JTJ) is written as H = Ho + V nT= i, where 

V^.=i=5^%cos(*a:i)®(|*><0| + |0)(*|). (3) 

i 

To gain intuition regarding the behavior of the system, let us define the operator V(x) as the optical potential 
operator V nT= i for which the position operators are replaced by definite positions x. In the N+l dimensional space of 
internal states for the ut = 1 manifold, the operator V(x) has two non-zero eigenvalues, ±Tigx(x) — ±Tig y / ^^cos^fcxi 
with corresponding eigenstates 

|D ± (x)) = -^^|0)±^£cosfea ;i |^. (4) 

We will refer to the |D_(x)) and |D+(x)) eigenstates of the potential matrix as the red and blue internal states, 
respectively, in reference to their energies being red- or blue-detuned from the empty cavity resonance. The remaining 
N — 1 eigenvalues of the ontica.l notentiai matrix axe null-valued. These corresnond to dark states havine no overlan 
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with the excited cavity internal state, |0), and which, therefore, cannot be excited by the cavity excitation interaction 
Hi. Note that (N c ) = 1/2(0) for all bright (dark) states, hence the cavity transmission spectrum is equivalent to the 
excitation spectrum in this treatment. We can now write the optical potential operator V nT= \ as 

V nT =i =9 J dx X (x) |x) (x| ® (|U+(x)) (D+(x)| - |U_(x)) (U-(x)|). (5) 

We also note that the initial state |<J//) can be written as a superposition of bright states, 

1*,) = -L ( |0 /(x ) ® D _(x)) + |^(x) ® £>+(x)» . (6) 

Our treatment allows us to recover easily results of the Tavis-Cummings model \2\\ in which a collection of fixed 
two-level atoms are coupled to a single-mode cavity with fixed, identical dipole coupling. Considering V( x o) with 
all atoms at the origin (xo = (0,0, . . .0)), we find a spectrum composed of delta-functions at ±g^/N (see Figure la) 
corresponding to the two bright states |D±(xo)). The clear dependence of the frequency of peak transmission on 
the integer number of atoms in the cavity provides the background for a basic, transmission-based atom-counting 
scheme. "Extrinsic" line-broadening, due to cavity decay and other losses, will smear out these sharp transmission 
peaks (see Figure lb), and will determine the maximum number of atoms that can be counted at the single- atom 
level by discriminating between the transmission spectra for N and TV + 1 atoms. For the remainder of the paper, we 
focus on intrinsic limitations to atom counting, i.e. those due to atomic localization and motion. 



III. METHOD OF MOMENTS 



To analyze the transmission characteristics of the atoms-cavity system in the presence of spatial dependence and 
atomic motion, we shall assume that the key features of the spatially-independent limit discussed above are maintained 
(Figure 2). Specifically, the transmission spectrum will still be described by two sidebands, one red-shifted and 
one blue-shifted from the empty cavity resonance by some frequency on the order of g. In determining the cavity 
transmission I(uj), we may thus divide the bright excited states {|^)} of the ut = 1 manifold into "red" — )} 
and "blue" states. From these "red" and "blue" states, we determine the transmission lineshapes I-(u>) and 

I+{oj) of the red and blue sidebands, respectively. 

The validity of this approach is made more exact by the following considerations. We have already obtained the 
locally-defined internal-state eigenbasis for the nr = 1 manifold as eigenstates of the operator V(x), namely the states 
\D± (x)) and the remaining N— 1 dark states. Let U (x) be the rotation operator which connects the uncoupled internal 
states {|0) , |1) , . . . \N)} to the eigenstates of V(x) at a particular set of coordinates x (the "coupled internal state 
basis"). Now, consider applying this local choice of "gauge" everywhere in the system. Since the dipole interaction 
operator V is diagonalized in the coupled internal state basis, it is convenient to examine the full Hamiltonian H 
in this basis. Defining the spatially-dependent rotation operator U = J dx |x) (x| U(x), we therefore consider the 
transformed Hamiltonian H' = UHW . 

Returning to Eq. the only portion of the Hamiltonian H which does not commute with the operator U is the 
kinetic energy. Considering the transformation of the momentum operator for atom i, 

Up&fi = Pl + -U^-U^ = Pi + A, (7) 

i aXi 

the transformed Hamiltonian H' can be expressed as H' = H a d + A_ff , where 



H - d = Y,{^® I+ 2 mLu2 ° X *® 1 ) + h 9X(x)(\D + )(D + \-\D_)(D_\) 1 (8) 
AH = ^-^(piAi + AiPi + AiAi). (9) 

i 

The operator H a d describes the behavior of atoms which adiabatically follow the coupled internal state basis while 
moving through the spatially-varying cavity field, and AH represents the kinetic energy associated with this local 
gauge definition. 
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Let us treat AH as a perturbation and expand the eigenvalues and eigenstates of H' as 



j? - p(°) 



I* 



(o) 



3,± 



(10) 
(11) 



We define projection operators onto the red and blue and dark internal states, n_,n+,ILj, respectively, with the 
explicit forms 



IE 



dx|x) (x|® |£> ± (x)>(£>±(x)|. 



(12) 



These projection operators commute with H a d- Hence the bright eigenstates of H a d, which are simultaneous eigen- 
states of II± and lid, can be written as 



D ± j = / dx0W(x) |x) ® |£>±(x)>. (13) 

We now assign an eigenstate, l^j), of i?' to the red or blue sideband if its zeroth order component belongs 

respectively to the \D_) or \D + ) manifold. We can therefore define the sideband transmission spectra I±(oj) as the 
separate contributions of red/blue sideband states to the total transmission spectra (see Eq. @): 



I ± (lo) cx K*j,±I*/)| 2 K u i,± - Wb ~ w). 
i 



(14) 



Determining the exact form of i±(w) is equivalent to solving for all the eigenvalues hu>j^± of the full Hamiltonian. 
This is a difficult problem, particularly as the number of atoms in the cavity increases. In practice, given the potential 
extrinsic line-broadening effects which may preclude the resolution of individual spectral lines, it may suffice to simply 
characterize main features of the transmission spectra. As we show below, general expressions for the various moments 
of the spectral line can be obtained readily as a perturbation expansion in AH. These moments allow one to assess 
the feasibility of precisely counting the number of atoms contained in the high-finesse cavity based on the transmission 
spectrum. 

In general, we evaluate averages (w+) weighted by the transmission spectral distributions I±(w). We make use of 
the straightforward identification (for notational clarity, shown here explicitly for the case of the blue sideband) 



% J did 1+ (id) id 

f(kdl+(id) 
E/V- '1'/ (ttj.+ l*/) 

Ei<*/| 1 */) 

E jt+ (n+ + n_) | 1 (n+ + n_)|* J 



Ej ( n + + n -) l*i,+> ( n + + n_) |*, 

where we have made use of the facts that 11+ + II _ +ILj = / and ILj \ tyj) — 0. To zeroth order, Eq. I|15|) becomes, 



(15) 
(16) 
(17) 



h(u + ) 



(0) 



2(* J |n + fr a<I ii + |*j) 



The first-order correction to this result is given by, 

= 2((* / |n+Affli+|* / ) + (* 7 |n_E,£ 



-4{<b I \n + H ad u + \^ I ) ((*/|n_Ei 



(0) 

(i) 



(i) 



(0) 



j-+ 



(0) 



n+|* i ) + <* J |n+ E,-4+ 



n+i^ + ^in+E,- 



(1) 



(0) 



(1) 



3,+ / \ 3,+ 

n_|*/ 



(18) 



n_ I*,; 



(19) 



To evaluate the sums over the first-order corrections to the eigenstates, 



(i) 



we approximate the energy denom- 



inator in the first-order perturbation correction as the difference between the average energies of the red and blue 
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AH 

*,|n_£|i#>)(*$iM*/> = EE 1 „ re) ' J ' +/ (*S n + !**> . ( 2 °) 



sidebands 

l,T,(°) \/,t,(°) AW ii 

J*L /,T,(0) 

1 



*/|n_Ai?n+ . (2i) 



Using this approximation, we can evaluate Eq. 115(1 to the first order in perturbation, yielding 

h{uj + ) = 2 (n + ini + ) + ^(o^^coa (4 (n + g ad Aim_) - 8 (n + H ad n + ) (n_AM + )), (22) 

where all expectation values are calculated over the initial state We can also calculate the second moment of 

the distribution using the same technique. To first-order, we obtain, 

h 2 (w%) = 2 (H + {Hl d + AHH ad + H ad AH) H + ) + ^-^-^ (4 <H + ff a 2 d AIffl_) 

-8(n + tf o yi + )<n_Afra + )). 

In order to evaluate these expressions, we must calculate expectation values of the form II ± H ad AH k U ± over the 
initial state To simplify matters, we note that we can act with the projection operators on the initial state \^i), 

which is equivalent to operating in the \D±) internal state basis. Since H ad is diagonal in the \D±) basis, and 4>i(x) 
is the A^-dimensional harmonic oscillator ground state, it is straightforward to obtain 

H ad |^(x)D±) = (E + %xW) |0i(x)U±) . (24) 
Using the definition in Eq. we find that the \D±) matrix elements of AH are given by the matrix 

^ U 2 



r/c 2 c(x) / 1 -1 



AH = —^® -1 1 )< ^ 



where we have defined, 



C(x) = -^ + l-E^^. (26) 



Combining Eqs. (|2*2*|) and with Eqs. and lt2*5)) . we obtain to first-order in AH, 

h ( U+ ) E = hg ( X ) + \^ (0 + {- <Cx) + (0 <X» , (27) 



(H) - (- + > 2 ) - *V (<X 2 ) - (xf) + ng?± (± «Cx 2 > + (C) <X 2 » - 2 (C) (X)) • (28) 

Here all expectation values are taken over the spatial state </>/(x). Although the function </>/(x) is simply the product of 
N harmonic oscillator ground states, the presence of various powers of x(x) and C( x ) m the above expectation values 
makes their analytic evaluation very difficult for arbitrary iV.To determine the dependence of these integrals on atom 
number N, one may expand the integrand as a Taylor series in x 2 , leading to approximate analytic solutions for the 
integral as a series in 1/N. After some tedious algebra, we find the average positions of the red- and blue-transmission 
sidebands to be 

hM-Eo = ±n 9 VN^(l-i^f)-^(^)+0(j ! ). (29) 

Here we quantify the relative length scales of the initial harmonic trap as compared to the optical interaction potential 
through the parameter e = exp (— /c 2 cr 2 ), which is related to the Lamb-Dicke parameter rj by V2rj — ka and a = 

Next, we obtain an expression for the width of the red and blue sidebands by evaluating the second moment of the 
sidebands. Expanding Eq. 1(28(1 as a series in 1/N, we obtain 

f? (H ) - <^ ± ) 2 ) = ±%V (1 - ef (1 + e) ± %^ - 7 =J— (1 - ef (3 + e) + (£) . (30) 
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To gain some physical insight into these results, we consider two important regimes: the tight and loose trap regimes. 
These different regimes are reflected in the corresponding values of the parameter e, which tends towards 1 in the 
extreme tight-trap limit and to in the extreme loose-trap limit. In the tight regime, the length scale of the trapping 
potential is much smaller than the wavelength of the light, i.e. ka -C 1. This is equivalent to the Lamb-Dicke regime 
and is applicable to current experiments for trapped ions in cavities 0, lldj , or for neutral atoms held in deep optical 
potentials [l^. I n the loose-trap regime, ka > 1 and atoms in the ground state of the harmonic oscillator potential 
are spread out over a distance comparable to the optical wavelength. As atoms in this regime sample broadly the 
cavity field, one expects, and indeed finds, a significant inhomogeneous broadening of the atoms-cavity resonance. 

In the extreme loose-trap limit (e — ► 0), we find 



(w±) - E /h = ±gJ- 1 - — + - — + O - , (31) 




<(A W± ) 2 ) = lg*±gl£--J£= + [i=). (32) 



In the loose-trap limit, the center of the red sideband is now located at gwN/2 instead of at gvN as we obtained for 
the spatially independent case. This difference is due to the spatial dependence of the standing mode; the atoms no 
longer always feel the full strength of the potential, but are sometimes located at nodes of the potential. We also see 
that the sidebands have an intrinsic width of w g/y/8. This width will play an important part in limiting our ability 
to count the number of atoms in the cavity in the limit of a loose trap. 

Considering the tight-trap limit, we expand in the small parameter ka and obtain 



(w±) - Eo/h = ±gVN (l - ifcV 2 j - \^k 2 a 2 + O (fcV) , (33) 

((A W± ) 2 ) = i 5 2 fcV±^^ifeV 4 + 0(fcV). (34) 
NV ' 1 8 y y 2m 2-JN 



In the limit ka — ► 0, the atoms are confined to the origin and we recover the Tavis-Cummings result discussed earlier, 
wherein the transmission sidebands are delta functions at ±g-*/N away from the empty cavity resonance. As the 
tightness of the trap decreases, the atoms begin to expericence the weaker regions of the optical potential and the 
centers of the sidebands move towards the origin. In addition, the sidebands develop an intrinsic variance which scales 
as k A a A . 

An important feature of both regimes is the intrinsic linewidth of both the red and blue sidebands (see Figure 2a) . 
This linewidth has a magnitude of approximately g^J (1 — e)/8 when the vacuum Rabi splitting is much larger than 
the atomic recoil energy, i.e., g 3> hk 2 /2m. It is unrelated to linewidth due to cavity decay or spontaneous emission 
which we have not addressed here and results purely from the spatial dependence of the atom-cavity coupling. Thus, 
it will provide an intrinsic limit to our ability to count N atoms, regardless of the quality of the cavity that is used. 
Our expression for the intrinsic linewidth also highlights an asymmetry between the red and blue sidebands. To 
first-order, increasing the atomic recoil energy reduces the linewidth of the red sideband but increases the linewidth of 
the blue sideband. Consequently probing the red sideband of the atoms-cavity system rather than the blue sideband 
would facilitate counting atoms. In addition, these results suggest that the ability to tune both the atomic recoil 
energy hk 2 /2m and the coupling strength g (this can be done, for instance, using CQED on Raman transitions) would 
be beneficial. We attribute the asymmetry between the sidebands to the different effective potentials seen by states 
within the red and blue sidebands. A detailed analysis of this aspect will be provided in a future publication. 



IV. CONCLUSIONS 



We have found that the transmission spectrum of the cavity containing N atoms trapped initially in the ground 
state of an harmonic potential will consist of distinct transmission sidebands which are red- and blue-detuned from 
the bare-cavity resonance, when the vacuum Rabi splitting dominates the atomic recoil energy. Analytic expressions 
for the first and second moments of the transmission sidebands were derived, and evaluated in the limits of tight and 
loose initial confinement. These expressions include terms containing the vacuum Rabi splitting hg and the recoil 
energy h 2 k 2 /2m. The former can be regarded as line shifts and broadenings obtained by quantifying inhomogeneous 
broadening under a local-density approximation, i.e. treating the initial atomic state as a statistical distribution of 
infinitely massive atoms. The latter quantifies residual effects of atomic motion, in essence quantifying effects of 
Doppler shifts and line broadenings. 
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These results can be applied to assess the potential for precisely counting the number of atoms trapped in a high- 
finesse optical cavity through measuring the transmission of probe light, analogous to the work of Hood et al. and 
Miinstermann et al. 0, for single atom detection. To set the limits of our counting capability, we assume that 
atoms are detected through measuring the position of the mean of the red sideband. In order to reliably distinguish 
between N and N + 1 atoms in the cavity, the difference between the means for N and N + 1 atoms must be greater 
than the width of our peaks, i.e., \(oj±(N)) — (ll>±{N + 1))| > Au± (see Figure 3). Let us consider that, in addition 
to the intrinsic broadening derived in this paper, there exists an extrinsic width k' due to the finite cavity finesse and 
other broadening mechanisms. Evaluated in the limit g 3> hk 2 /2m and assuming large N, 



W-(-Y)} -U \ "A/^r- (35) 



We thus obtain an atom counting limit of 



ip + I(l-e)2(l + e) 



(36) 



where we have assumed that the intrinsic and extrinsic widths add in quadrature. This atom counting limit ranges 
from N max = g 2 /An 12 in the tight-trap limit, to N max = 1/(1/2 + 8k' 2 / g 2 ) in the loose trap limit. Figure 4 shows 
N max as a function of e for various values of k. In general, atom counting will be limited by extrinsic linewidth when 
16k' 2 > g 2 (l — e) 2 and by intrinsic linewidth when 16k' 2 < g 2 (l — e) 2 . 

These results demonstrate that atom counting using the transmission spectrum is best accomplished within the 
tight-trap limit. Certainly, in the loose-trap limit, atom counting will be rendered difficult as the intrinsic linewidth of 
the sidebands is increased. However, several questions regarding the feasibility of atom counting experiments remain. 
First, although atom counting by a straightforward measurement of the intensity of the transmitted light may be 
difficult, it is possible that the phase of the transmitted light may be less affected by motional effects 0|. Dynamical 
measurements (possibly using quantum feedback techniques) might also yield higher counting limits. Second, atomic 
cooling techniques could be used in the loose-trap limit to cool the atoms into the wells of the optical potential, 
thereby decreasing the observed linewidth [2^, |2|| ^ • Finally, the state-dependence of spontaneous emission has 
not yet been taken into account. Although the loose-trap regime leads to an intrinsic linewidth which limits atom 
counting, it may also suppress the extrinsic linewidth as a result of contributions from superluminescence. On the 
other hand, in the Lamb-Dicke limit, the atoms are all highly localized, which could lead to enhanced spontaneous 
emission due to cooperative effects. Future work will investigate alternative methods of atom counting and will explore 
complementary techniques of reducing the intrinsic linewidth in atom-cavity transmission spectra. 
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FIG. 2: a. Intrinsic transmission spectrum of atoms-cavity system including spatial dependence of potential and atomic motion, 
b. Corresponding transmission spectrum including cavity decay. 
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FIG. 3: Plot of (uj-) as a function of the trap tightness e = exp(— k' 2 a 2 ) for N = 8 and N = 9 and small ratio of atomic 
recoil energy to vacuum Rabi splitting, Tik 2 /2mg = 0.01. The shaded regions indicate the instrinsic width of the red sideband, 
±a/ ((Au-) 2 } jl. In the tight-trap limit, TV = 8 and N = 9 can be distinguished. In the loose-trap limit, the intrinsic width of 
the spectra render determination of atom number difficult. 
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FIG. 4: Maximum limit N max on atom counting as a function of trap tightness e — exp (— k 2 a 2 ) for several values of the decay 
parameter k. e — > corresponds to the loose-trap limit while e — > 1 corresponds to the tight-trap limit. Notice that for the 
infinitely tight trap, atom counting is limited only by k. 



